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ABSTRACT 

We consider the effect of lensing by galaxy clusters on the angular 
distribution of submillimeter wavelength objects. While lensing does not change 
the total flux and number counts of submillimeter sources, it can affect the 
number counts and fluxes of flux-limited samples. Therefore imposing a flux 
cut on point sources not only reduces the overall Poisson noise, but imprints 
the correlations between lensing clusters on the unresolved flux distribution. 
Using a simple model, we quantify the lensing anisotropy induced in flux-limited 
samples and compare this to Poisson noise. We find that while the level of 
induced anisotropies on the scale of the cluster angular correlation length is 
comparable to Poisson noise for a slowly evolving cluster model, it is negligible 
for more realistic models of cluster evolution. Thus the removal of point sources 
is not expected to induce measurable structure in the microwave or far-infrared 
backgrounds. 



Subject headings: gravitational lensing - galaxies: clusters: general - infrared: 
galaxies - cosmic microwave background 



- 2 - 



1. Introduction 

Due to the steep slope of the Rayleigh- Jeans portion of their emission spectra, dusty 
galaxies have submillimeter luminosities that are greatly enhanced as a function of redshift. 
This results in large negative submillimeter K-corrections, causing the number of galaxies as 
a function of observed flux to be roughly constant between redshifts 1 and 10 (Franceschini 
et al. 1991; Blain & Longair 1993). Thus observations that can resolve submillimeter 
sources at z — 1 should be able to detect the full history of these objects, and one would 
expect submillimeter source counts to have a much steeper slope as a function of limiting 
flux than source counts in other wavebands. 

While this strong dependence on limiting flux had long been predicted (Blain & 
Longair 1993; Blain & Longair 1996), the overall surface density of these objects remained 
uncertain to within a factor of a thousand until recently. Using the the new Submillimeter 
Common User Bolometer Array (SCUBA) (Holland et al. 1999) on the James Clerk Maxwell 
Telescope, several groups (Barger et al. 1998; Eales et al. 1993; Hughes et al. 1998; Smail, 
Ivison, & Blain 1997) have detected objects at 850 /iin with fluxes of a few mJy Meanwhile, 
parallel investigations utilising the FIRAS instrument on COBE (Puget et al. 1996; Fixsen 
et al. 1998; Hauser et al. 1998), DIRBE combined with IRAS (Schlegel, Finkbeiner and 
Davis 1998) and TeV gamma ray constraints via pair production against the IR background 
(Biller et al. 1998) have begun mapping out the extra galactic background from ~ 5 to 
1000 fim wavelengths. Thus the past few years have not only given us our first look at the 
high-luminosity submillimeter universe, but allowed us to place strong constraints on the 
unresolved distribution as well. This has allowed several authors to investigate the effects 
of confusion noise from unresolved sources both on galaxy surveys and measurements of the 
cosmic microwave background (CMB) (Blain, Ivison, & Smail 1998; Blain et al. 1998; Scott 
& White 1998). 

The steepness of submillimeter source counts also makes them highly sensitive to the 
effects of gravitational lensing. As a gravitational lens amplifies the flux of objects behind 
it while at the same time reducing the effective volume that is being seen, only objects 
with sufficiently steep source counts display a statistical increase in bright-end counts when 
lensed. Thus one would expect the lensed fraction of bright submillimeter-wave sources to 
be much larger that the fraction in the optical and radio wavebands (Blain 1996). This 
sensitivity has motivated several SCUBA searches for objects behind galaxy clusters, using 
lensing to boost resolved source counts. It has also led to investigations of the effect of 
lensing by both galaxies and galaxy clusters on overall number counts and the prospect 
of measuring the degree of lensing from all-sky microwave experiments (Blain 1996; Blain 
1997; Blain 1998). 



3 



In addition to changing the number counts of bright objects, lensing induces structure 
in the unresolved source distribution. As the clumping of matter cannot affect the total 
flux or number of sources, lensing skews the flux distribution by making certain sources 
brighter while at the same time preserving the total flux and number of objects. Thus 
imposing a flux cut causes a larger flux decrement in the areas behind gravitational lenses. 
As clusters of galaxies are known to have large angular correlations (Dalton et al. 1984; 
Tadros, Efstathiou, & Dalton 1998), the removal of resolved point sources will imprint a 
correlated decrement on the unresolved distribution of submillimeter sources. 

As future experiments such as the Planck satellite (Bersanelli et al. 1996) will be 
making sensitive all-sky measurements at microwave and submillimeter wavelengths and will 
use flux-limited samples to investigate the structure of background emission, it is important 
that this effect be quantified. While many elements contributing to this structure remain 
uncertain, current cluster and submillimeter observations have improved our knowledge to 
the point where this quantification is now possible. 

In this work we develop a simple model of lensing by clusters, appropriate to the 
degree at which cluster evolution is understood. We then combine our lensing model with 
a semi-empirical model of the submillimeter universe (Tan, Silk & Balland 1999, hereafter 
TSB), based on recent observations. This leads to a quantification of induced structure by 
lensing that is useful in determining the extent to which this effect needs to be considered 
in future observations. 

The structure of this work is as follows: In Section 2 we discuss various models of 
the lensing of high redshift objects by clusters. In Section 3 we examine the number and 
spectral evolution of the submillimeter sources themselves and apply our lensing models to 
predict the number counts, fluxes, and clustering properties of lensed sources. In Section 
4 we compare these predictions to Poisson noise and discuss implications for background 
experiments. Conclusions are given in Section 5. 



The effect of gravitational lensing on the number densities of distant sources can be 
described by F(A, z)dA, the probability that a source at redshift z will have its flux density 
amplified by a factor between A and A + dA. For any distribution of objects 



2. Models of Gravitational Lensing 



2.1. Lensing by Clusters 




(1) 
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the first equality being a statement of conservation of probability, and the second following 
from flux conservation. Note that while naively one would expect lensing to cause only 
amplification, the requirement that the mean flux density be unaffected by the clumping 
of matter in the universe results in a net deamplification of unlensed objects (Weinberg 
1976; Dyer & Roeder 1973). Thus for all redshifts F(A,z) is a skewed distribution with a 
maximum slightly less than one and a long tail corresponding to strong lensing events. 

By integrating over the comoving mass distribution of clusters, n(M,z), F(A,z) can 
be determined as 

where dr is a comoving radial distance element, and b is the proper impact parameter for a 
source at redshift z and a lens mass M located at a redshift z' (Peacock 1982). 

The proper impact parameter is the angular size distance to the lens times the angle 
of deflection by that lens, and is dependent on the geometry of the lensing distribution. 
Fortunately, the most simple assumption, the p oc r~ 2 isothermal sphere distribution, has 
been shown both theoretically and observationally to be a good fit for intermediate radii 
between 0.1 and 1.5 h^ 1 Mpc, where h is the Hubble constant in units of 100 kms" 1 Mpc -1 
(Tyson & Fischer 1995; Navarro, Frenck & White 1996). For large magnifications this 
model gives (Peacock 1986) 

were a v is the velocity dispersion of the lens, D(0,z') is the angular size distance to the 
lens, D(0, z) is the angular size distance to the lensed object, and D(z',z) is the angular 
size distance from the lensed object to the lens. For clusters, a v = 860 m^km/sec, where 
m is the mass of the cluster within 1.5Mpc/i _1 in units of 6 x 10 14 M Q /i -1 . Strictly 
speaking, the angular size distances in Eq. |3| are the Dyer-Roeder distances which account 
for the inhomogeneity of the lensing objects. For our purposes here, however, we use the 
homogeneous distances (Young et al. 1980). 

Bahcall, Fan, & Cen (1997) have used the CNOC/EMSS sample of high-redshift 
clusters to determine the observed evolution of cluster abundances. For massive 
clusters with m > 1.05, the evolution in redshift is given by \ogn(z) oc — z with 
log(n(z = 0)/n(z = 0.5)) = 0.7 ± 0.3(1<t) (Fan, Bahcall & Cen 1997, hereafter FBC). 
Carlberg et al. (1996) give similar results. One could envision modelling the evolution of 
clusters from a Press-Schechter (Press & Schechter 1974) viewpoint, counting the number 
of peaks above a certain threshold as a function of redshift and cosmology. The problem 
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with this approach, however, is that while the evolution of the number density is only 
weakly dependent on the cosmological parameters, it is strongly dependent on the unknown 
amplitude of mass fluctuations on the 8 h~ x Mpc scale, <7g (FBC). As og is uncertain 
from observations, a simple fit to data is the best we can hope for even if we fix the other 
cosmological parameters. 

For our purposes then, we adopt a simple model for the number and redshift evolution 
of clusters that retains the Press-Schechter form but is fit to the observed number density 
and evolution of clusters 

^L = nom -V^ e -i^Hz a+ P)), (4) 
dm 

Here the m 1//2 dependence in the exponential and overall m~ 7 / 4 dependence corresponds 
to a mass density power spectrum P(k) oc fc~ L5 around 8 h^ 1 Mpc, where cold dark 
matter models give P(k) oc k~ 1A (Bardeen et al. 1986) We fix these parameters to the 
n(> M) relation for clusters as calculated in Carlberg et al. (1997), which used the 
Canadian Network for Observational Cosmology (CNOC) (Yee, Ellingson, & Carlberg 
1996) cluster sample drawn from the Extended Medium Sensitivity Survey (EMSS) (Henry 
et al. 1992) as well as the ESO Cluster Survey (Mazure et al. 1996) and other cluster 
data sets (Henry & Arnaud 1991, hereafter HA; Eke, Cole, & Frenck 1996, hereafter 
ECF). This comparison yields (3 = 10 and n x e _/3 = 10 _6 /i~ 3 Mpc _3 . For this set of 
parameters the log(n(z = 0)/n(z = 0.5)) comparison of FBC quoted above corresponds to 
a = 2.9 ± l.l(lcr). Note, however, that cluster evolution fits have been carried out by other 
authors and there is some disagreement between results (Blanchard & Bartlett 1998; Sadat, 
Blanchard, & Oukbir 1998). In Fig. [1] we compare the data compiled in Carlberg et al. 
(1997). with the "best-fit" FBC, a = 2.9 evolution model and a slowly evolving, a = 0.7 
model that we consider in further detail below. 

Modelling clusters in this manner, we find that for large amplifications 

F(z, A) = a(z)A~ 3 where 

a(z) = 0.007 r dzl ^SWWS(n(z-*)r p» W / V ° e -<">''"<-+">>>, (5) 

JO dz' b{T](Z)) A Jm min 

where rj is a dimensionless comoving distance (77 = rH/c), S(rj) = rj in a flat universe, 
and S(r]) = smh(rj\/l — Qm — ^a)/\/1 — — in an open universe. Note that this 
expression is comparable to the model considered by (Pei 1995), which he derived from 
observations of the empirical distribution of dense groups and rich clusters (Zabludoff 
et al. 1993), although our approach is quite different. Throughout this paper we fix 
M min = 3 x 1O 14 M and M max = 1 x 1O 15 M . 



In Fig. we plot this function for three different cosmologies in the upper panel, and 
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Fig. 1. — The n(> M) relation as calculated fit to the data compiled by Carlberg et al. 
(1997), with redshift ranges indicated and acronyms are as described in Sec 2.1. The solid 
lines are z — 0, the dotted lines are z = 0.5, and the dashed lines are z = 1.0. The 
upper panel shows the FBC (a = 2.9) model and the lower panel shows the slowly-evolving 
(a = 0.7) model. M is the mass of the cluster within 1.5k" 1 Mpc. 
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Fig. 2.— Upper panel: Plots of a(z), a measure of the probability of strong lensing as 
function of redshift, for three different cosmological models: a flat, Qm = 1 universe (solid), 
an open model with Qm = 0.3 (dotted), and a flat cosmological constant model with Qu = 0.3 
(dashed). In all cases a = 2.9. Lower panel: Plots of a(z) for Einstein-de Sitter cosmologies 
with four different models for the evolution of clusters: a = 4.0 (dotted), a = 2.9 (solid), 
a = 1.8 (short-dashed) and a = 0.7 (long-dashed). 
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for four different cluster evolution models in the lower panel. Notice that at moderate 
redshifts the cosmological constant model gives slightly more than twice the lensing of the 
Einstein-de Sitter model, while the slow-evolution model of lensing is four times higher 
than the FBC evolutionary model. Thus while the degree of lensing is quite sensitive to 
cosmological parameters (Carroll & Press 1992), lensing by galaxy clusters is even more 
sensitive to the evolution of the cluster population itself. 



2.2. Clustering of Clusters 

The effect of gravitational lensing on the clustering of a distribution of sources depends 
on the probability that two sources at redshift z\ and z 2 separated by an angle 9 will be 
lensed by factors between A\ and A\ + dAi and A 2 and A 2 + dA 2 : F(A\, z%, A 2 , z 2 , 9)dAidA 2 . 
Similar probability and flux conservation relations to Eq. p] hold for this distribution: 

/•OO /"OO /"OO /"OO 

/ dA 1 dA 2 F{A 1 ,z 1 ,A 2 ,z 2 ,9) = l = A l dA l / A 2 dA 2 F(A u z u A 2 , z 2 , 9). (6) 

J J J J 

In this case F(A±, z±, A 2 , z 2 , 9) can be determined as 

F(A 1 ,z 1 ,A 2 ,z 2 ,9) = / — dz [{1+z 2 x 

dA\ dA 2 Jo dz[ 

^jdz' 2 (l + ztfnbiA^nbiArfn 2 ^, z 2 , 9), (7) 

U / Z 2 

where n 2 (zi, z 2 , 9) is the product of the number density of clusters at redshifts Z\ and z 2 
separated by an angle 9. 

Strictly speaking 

n 2 {z u z 2 ,9) =n(z 1 )n(z 2 )[l + £(r(z 1 ,z 2 ,9))}, (8) 

where ^(r) is the cluster spatial correlation function, and r(zi, z 2 , 9) is the comoving distance 
between two clusters as a function of their redshifts and separation. There are however 
a number of difficulties in implementing such an approach, as not only is r dependent on 
unknown cosmological parameters, but the spatial correlation function of clusters is subject 
to a number of observational uncertainties. While several authors have fit the optical 
cluster-cluster correlations with a power law of the form 7] = (r/r ) -7 with —2.2 < 7 < —1.8 
(Peacock & West 1992; Dalton et al. 1994), the cluster-cluster correlation length ro has 
remained controversial. Here values are much more uncertain, ranging from 15 h^ 1 Mpc 
(Nichol et al. 1992; Dalton et al. 1994) to 25 h' 1 Mpc (Peacock & West 1992; Abadi, 
Lambas, & Muriel 1998), with various dependences on cluster mass and mean intercluster 
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separation (Bahcall & Soneira 1983; Bahcall & Cen 1992). To compound this problem, the 
clustering of clusters should be more severe at higher redshifts due to the fact that the first 
density peaks to collapse were more rare than those collapsing recently. 

Due to these uncertainties, we avoid implementing the spatial correlation function and 
approximate the number density of objects as 



n 2 (z l , z 2 , 9) « n(z l )n(z 2 )[l + w 



(9) 



with w(6) = (6>/6> ) 1+7 , with O = 1.1° observed in the APM survey (Dalton et al. 1994). 
Note however, smaller values of 9 may be more appropriate for lensing due to the larger 
redshift range at which clusters have an observable effect. Thus our choice for the angular 
correlation function of lensing clusters should be considered an upper limit. 

As the form of the peak of F(A, z) at amplifications slightly less than unity has little 
effect on the results of lensing, Blain (1996), following Vietri & Ostriker (1982), examined a 
simple model for F: 



F = H(z)5(A - A (z)) + A-*a(z)6(A - A min )9(A 



A), 



(10) 



where H(z) and A(z) are fixed by the conservation relations Eq. [1|, and 9 is the Heaviside 
step function, which is used to exclude low magnifications at which we would expect Eq. ^] 
to break down and high magnifications at which the isothermal sphere approximation is 
invalid. Working to first order in a(z), this gives 



H(z) 



1 —a(z) 



A' 2 - A~ 2 

^min ^max 



A {z) 



1 — a(z) 



9 A ■ — 1 

?A 2 ■ 



94 ■ — 1 
^-^min x 

2A 2 

max 



Physically, this corresponds to modelling each cluster as an isothermal sphere in which we 
have truncated the lens profile at the angular size distance corresponding to an amplitude 
of A min , smoothed the distribution of matter within a core radius corresponding to A max , 
and imposed a net demagnification of objects relative to their amplitudes in a homogeneous 
universe. 

In order to study the clustering properties of lensed objects, we examine a similar 
model for the two point amplification function 



F(A 1 ,z 1 ,A 2 ,z 2 , 



= H( Zl , z 2 , 6)5{A 1 - Aofa, z 2 , 0))5{A 2 - A ( Zl , z 2 , 9)) 
+ AfAf9(A l - A min )9(A 2 - A min ) x 

0(A max - A 1 )9(A miiX - A 2 )a{z 1 )a{z 2 )[l + w 



(12) 
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where H and A are fixed by conservation of probability at all values of redshift and angular 
separation. Imposing these conditions gives us 



H(zl,z2,9) = l-a(z 1 )a(z 2 )[l + w(9)} 
A Q (zl,z2,e) = l-a(z 1 )a(z 2 ) l+ ^ 9) 



9/42 o/l 2 
min max 



A ■ A 

-^min ^iiiax 



1 



1 



2AL„ 2AI 



(13) 



Fixing the typical unlensed size of a submillimeter source to be 25/i kpc and the 
redshifts of a typical cluster and submillimeter source to be 0.5 and 3 respectively, gives 
a value of A max = 70 in a flat universe (Qm = 1), with higher values for open cases. 
Throughout this work we fix A max = 50, as a conservative limit. Somewhat more arbitrarily, 
we have set a minimum threshold of A m[n = 2, because it is a round number above which 
more sophisticated lensing models typically display F(A, z) oc A~ 3 profiles (Peacock 1982; 
Pei 1995). We have examined the effect of these cutoffs on our model of cluster lensing and 
find that changing them introduces less than 20% uncertainty in our results for reasonable 
value of A max and A min and relevant threshold fluxes. 



3. Submillimeter Sources 

3.1. A Semi-Empirical Model 

We use the simple model of TSB for the evolution of sources in the infra-red and 
submillimeter. The model evolves the properties of the local infra-red and submillimeter 
sources backwards in time, accounting for dust, gas, and spectral evolution. It agrees well 
with observations and constraints of the extra-galactic background from ~ 5 to 1000 /im 
(Puget et al. 1996; Fixsen et al. 1998; Hauser et al. 1998; Biller et al. 1998), as well as with 
number count surveys from ISO (Puget et al. 1999; Aussel et al. 1998; Altieri et al. 1998) 
and SCUBA (Smail et al. 1998). Two types of sources, disk galaxies and starbursts, have 
significant fluxes in the submillimeter regime due to the reprocessing of stellar radiation by 
dust. 

Disk submillimeter luminosity evolution is derived from observations of the star 
formation rate and metallicity histories of the Milky Way. A Schmidt Law relation 
(Kennicutt 1998) between star formation rate and gas density is used to derive an optical 
depth history. The contributions from both young and old stars to dust heating are 
modelled. The number evolution of disks is obtained from models of collision-induced 
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galaxy formation (Balland, Silk, & Schaeffer 1998) which attributes the formation of 
different morphological types to galaxy-galaxy interactions. 

We model starburst evolution assuming they result from the mergers of gas-rich 
disks. The typical luminosity increases with the gas fraction, but decreases with the mean 
mass of the galaxies. The starburst number density is proportional to the instantaneous 
galaxy-galaxy collision rate. 

These evolutionary models are applied to the present day far infra-red luminosity 
function, as observed by IRAS from 40 to 120 /im (Saunders et al. 1990), assuming an 
Einstein-de Sitter cosmology with h = 0.5. Given the uncertainties in the lensing model 
due to the evolution of clusters, we have not extended the submillimeter model to other 
cosmologies, but rather we restrict our attention to this cosmology. As lensing is more 
severe in open models and models with nonzero Q\, these are expected to yield results 
above the a = 2.9, = 1 model, although we still consider the a = 0.7, Qm = 1 model we 
examine below to be an overestimate of lensing even compared to Qm < 1 cosmologies. 

For simplicity, a cut is made in luminosity at L cnt ~ 6 x 10 10 h~ 2 L Q , below which 
sources are treated as disks, and above which they are treated as starbursts (Sanders & 
Mirabel 1996). The results show the most sensitivity to this parameter, and to illustrate 
model uncertainties, we investigate the effect of varying L cut in Sections 3.2 and 4. 

The luminosity function is also divided empirically into different spectral classes 
(Malkan & Stecker 1998), and these spectra are analytically extended beyond 400 /im into 
the submillimeter regime, assuming a dust emissivity dependence of v 13 . Here we fix (3 = 1.5 
(Franceschini, Andreani, & Danese 1998; Roche & Chandler 1993), but our results are quite 
insensitive to this value as most of the sources detected are at z > 1. 



3.2. Lensed Submillimeter Sources 



Applying Eq. [LT] and [II] to our model for number counts of submillimeter objects as a 
function of flux and redshift, we are able to calculate the perturbation to the number of 
sources per unit flux due to lensing by clusters. If we take rf ^° tal 
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(14) 
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where ^ = ^ a{z)-§^dz. Notice that both / °° dSA% and / °° SdSA% are zero, as 
the lensing of objects does not change the total number of observed objects or the total 
integrated flux from those objects. Suppose, however, we consider only objects above a 
certain threshold flux S cnt . fg° ut dSAjj^f and SdSAj^ need not be zero, as both the 
total number of objects above a given threshold and the total contribution to the flux by 
these objects can vary according to the degree of lensing. Furthermore, the total change 
above a threshold must be equal to the change below that threshold. Thus 



9 dSA lS 

Ocut LiLj 



Scut dN r°° dN /^cut dN 

dSA^A, and / SdSA^— = - I SdSA 

Scut do 



10 dS J s^t dS Jo dS 

for all values of S cnt , and we are free to consider these quantities interchangeably. 

d 2 N 2 



(15) 



By applying Eq. [T2| and [13], the perturbation to dSidS2 due to lensing can be similarly 
calculated. Here we find 
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d 2 N 2 
dSidSi 
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(16) 



We now restrict our attention to flat cosmologies with h = 0.5. Given the uncertainties 
in our knowledge of clusters and submillimeter sources, we expect this model to be 
sufficiently representative to examine the importance of structure induced by lensing. Due 
to the large volume and degree of lensing in open and cosmological constant models, we 
expect a = 2.9, < 1 models to fall somewhere between our Einstein-de Sitter a = 2.9 
and slowly evolving cluster models. 

In Fig. ^ we plot the total number of objects above a given threshold flux as a function 
of S cut as well as the change due to lensing by clusters at 850 /mi, roughly corresponding 
to the SCUBA observations and one of the high-frequency channels on the Planck satellite. 
The model shown in this figure is not fine tuned to match the sub-mm observations, 
but rather is one example of a range of possible models. While the number of strongly 
lensed objects is equal to AiV(> 5 cut ) for large values S CVLt , deamplification causes the net 
perturbation to iV(> S cnt ) to be an order of magnitude smaller than the number of strongly 
lensed objects at lower S cut values. 

As noted by Peacock (1982), gravitational lensing is most important when ^ is steeper 
than S~ 3 and is able to overcome the A~ 3 dependence of the high-magnification tail of 
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Fig. 3.— Number density of resolved objects as function of flux threshold at 850 /zm. 
The solid line is the total number density of objects, N(> S cnt ), from TSB. Note that his 
particular model is not fine tuned to match the sub-mm observations, but rather is one 
example of a range of possible models. The lower curves show the perturbation due to 
lensing. For the Qm — 1, a — 0.7 model the short-dashed line shows the change in number 
densities due to strong lensing by clusters, AiV(> S CVLt ), and the long-dashed line is the 
number of strongly lensed objects (the contribution due to the last term of Eq. 14). The 
dotted curve shows AiV(> S cnt ) for the FBC Einstein-de Sitter model with a = 2.9. The 
points are the SCUBA data as tabulated by Smail et al. (1998). All errorbars are la. 
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the lensing distribution. For more shallow distributions, however, the degree of lensing is 
smaller by several orders of magnitude, and can even result in a net decrement to N(> S cut ) 
for sufficiently shallow slopes (Broadhurst, Taylor, & Peacock 1995). Thus the difference 
between the number of strongly lensed objects and AiV(> S cu t) becomes important at 
values roughly corresponding to the turnover in the number counts of submillimeter objects. 
In the a = 2.9 and slowly evolving models AN(> S cnt ) is almost of the same form, but is 
smaller for the slowly evolving model by a factor ~ 3. 

In Fig. | we compare flux per unit angle due to unresolved objects in lensed and 
unlensed distributions as a function of flux cut at 850 /zm. As objects behind gravitational 
lenses are fewer and more luminous, imposing an overall flux cut causes a stronger decrement 
in the lensed distribution of objects than in the unlensed case. Again, this flux decrement 
in a slowly evolving model is of the same form as in the FBC model, but scaled down by a 
factor ~ 3. Note that while the flux of lensed objects is only ~ 0.1% of the total flux, the 
correlations induced by clustering can make this contribution important at large angular 
scales. 

This discrepancy in fluxes also raises the question of the clustering of the objects 
themselves. While this is likely to be weak and at very small angles, Scott & White (1998) 
have demonstrated that this clustering may be dominant over Poisson noise at angles ~ 2', 
although their model is only for illustrative purposes. As this contribution goes down as 
a function of flux-cut, and is important only at small angles, for our purposes it is best 
thought of as an uncertainty in the small-scale noise and will turn out not to affect our 
results. 

Our models are quite sensitive to the threshold between starburst and disk galaxies. 
Raising the value of L cnt suppresses the high-;? luminous sources, shifting the turnover 
point in the ^ curve. As the lensed flux is extremely sensitive to this slope, the 
number counts and flux decrements in these model differ dramatically. In Fig. |5| we 
consider two submillimeter models that overpredict and underpredict this slope. The 
steep (L cut = 3 x 10 10 h~ 2 L Q ) model displays a large peak in the number counts and flux 
decrements that is absent in the more shallow (L cnt = 12 x 10 10 h~ 2 L Q ) model. While these 
are extreme models, chosen to make a point, ignorance of this slope is nevertheless one of 
the largest uncertainties in our analysis. 
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Fig. 4. — Flux per unit steradian in unresolved objects as a function of flux threshold at 850 
/im. The solid line is unlensed flux, the dashed line is the decrement in flux due to lensing 
by clusters in the slowly evolving model (Qm — l,cx — 0.7) and the dotted line is the flux 
decrement in the Qm — l,a — 2.9, model. 
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Fig. 5. — The effect of uncertainties in the submillimeter model on lensing. In the upper 
panel, the upper curves show the total number of objects without lensing by clusters, 
N(> S cut ), for two different models while the lower curves show the change in this number 
density due to lensing, AN(> S cnt ). Here L cut = 3 x 10 10 h~ 2 L & is represented by the 
solid lines, and the dotted lines represent L cnt = 12 x 10 w h~ 2 L Q . In the lower panel, the 
upper curves show total unlensed fluxes and the lower curves show the flux decrement due 
to lensing, with the lines corresponding to models as in the upper panel. 
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4. Background Anisotropies 

Contributions of point sources to background anisotropies can arise in two different 
ways: Poisson noise due to the intrinsic graininess of the distribution, and structured noise 
due to the placement of the sources on the sky. The first source of noise is given by 

/•Scut „dNrr± r s °<* „ / dN dN\ 

where we use the notation C\ to denote the angular power spectrum in units of intensity 
rather than the usual decomposition in terms of the temperature fluctuations. C\ is related 
to the usual Cg by 

Ct = C{ (tcmb^) 2 = C\ ^2.71 x 10 10 mJy sr" 1 ^^) " , (18) 

where x = hu/k B T CM B = Z//56.84 GHz and 1 Jy = 1(T 26 Win" 2 ^" 1 (Scott & White 1998). 
Note that here we follow the conventional notation for microwave background fluctuations, 
expressing our results in terms of C^s, given by the Legendre transform of the correlation 
function C{6). 

As lensing skews the distribution by gathering the flux into fewer, brighter objects, 
it does slightly increase the Poisson noise of the observations. This contribution to the 
anisotropies is a single-point process, however, with no intrinsic scale, and is completely 
independent of the placement of the lensing objects on the sky. The level of Poisson noise 
in our model and the increment due to lensing are shown in Fig. || Notice that removing 
point sources above a certain flux only decreases the Poisson noise, as it reduces the total 
power in the grainy distribution. 

The second contribution of point sources to anisotropies, that of inhomogeneities in 
the mean flux as a function of position is given by 



/•Scut /"Scut 

C{6y = / SidS! / S 2 dS. 
Jo Jo 



A dN 2 (S 1 ,S 2 ,9) 
dSidS 2 



(19) 



where A ^j^f^' 8 ^ ^ s as S^ ven m ^q- |16|- Note that at this point we are ignoring intrinsic 
correlations of the submillimeter objects themselves. While this is a concern at high I 
values > 500 corresponding to the galaxy-galaxy correlation scale (Scott and White 1998), 
it introduces little uncertainty into the analysis at is corresponding to cluster correlations. 

If one imposes no upper limit on the flux of the point sources, then the lensing cannot 
add structure to the angular correlations of point sources. This is due to the conservation 
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of total flux behind the lenses, Eq. 15, which insures that statistically, the average flux 
from a lensed pixel will be equal to that of unlensed areas of the sky. Unlike Poisson noise, 
imposing a finite flux cut by removing point sources actually increases this contribution to 
anisotropies as it skews the mean flux from lensed pixels to lower values. This introduces a 
flux anisotropy where there previously was none, adding power on the angular scale given 
by the clustering of the lenses. 



We can calculate the magnitude of this effect from the Legendre transform of Eq. [T! 

2" 



c! 




(IS2S2 



dS 2 



Ai dS 2 



(S a / A) (20) 



A\ dS x 

where w £ « 2-kOq' 1 This expression is only at £ values below the angular size 

of a cluster at a typical redshift ~ 0.5, and at £ values above which a power law fit to 
the angular correlation function of galaxies is appropriate (5° > 9 > V or 10 < £ < 1000) 
(Tadros, Efstathiou, & Dalton 1998). 

In Fig. [7] we plot this contribution to the anisotropy of unresolved sources along with 
the Poisson noise anisotropy at £ = 1/(1.1°) = 27. As Cg is proportional to (^p-) 2 , the 
slowly evolving lensing model gives a contribution ~ 25 greater than the a = 2.9 model. 
Note also that this comparison is only valid at this particular £ value as the Poisson noise 
and lensing anisotropy have different scalings as a function of £, with £ 2 Ci of the lensing 
contribution going as £~ 7-1 , and £ 2 C% of Poisson noise going as £ 2 . Thus, at higher £ 
values, Poisson noise will be more dominant and will be even more enhanced at very high £ 
values by galaxy-galaxy correlations, while at lower values the lensing anisotropy will play 
a larger role. For reference, the expected microwave background contribution is given by 
£ 2 C £ /2n « 10- 10 for £ < 1000. 

In Fig. |H we plot the same quantities, but now for the L cut = 3 x 10 10 h~ 2 L Q 
submillimeter model. While this model slightly overproduces the background flux and most 
likely overpredicts both lensing and Poisson anisotropies, it nevertheless serves to illustrate 
the sensitivity of our results to the slope of the number counts near the flux-cut threshold. 
While Poisson noise is boosted only by a factor of two in this model, the lensing anisotropy 
is increased by more than an order of magnitude by a change in this slope. Nevertheless, 
even in this extreme model, the FBC lensing anisotropy falls short of the Poisson noise. 



Finally, in Fig. we compare Poisson noise and noise induced by lensing at 550 fim 
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Fig. 6. — Contribution of confusion noise to measurements of Ci/2tt as a function of flux cut 
The solid line is unlensed flux while the dotted line is the increment in this flux due to lensing 
by clusters. The dashed line is the Poisson noise as calculated from a simple double-power 
law fit by Scott & White (1998). 




Fig. 7.— Contribution of lensing-induced structure to measurements of l 2 Ce/2n as a 
function of flux cut at £ = l/9 , the correlation angle. The solid line is the Poisson noise 
of the unlensed distribution assuming 9 = 1.1°, the dashed line is the anisotropy due to 
lensing in the slowly evolving cluster model, and the dotted line is the anisotropy in the FBC 
evolution model. 
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Fig. 8. — Contribution of lensing-induced structure to measurements of £ 2 Ce/2Ti at £ = 1/0 . 
as a function of flux cut, for a model (Lcut = 3 x 10 10 h~ 2 L Q ) with a steeper source count 
slope near the flux cut The curves are as in Fig. [7|. 
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Fig. 9. — Contribution of lensing-induced structure to measurements of £ 2 Ce/2ii at £ = 1/6q 
as a function of flux cut at two different wavelengths. The upper panel corresponds to the 
550 /im, 545 GHz Planck channel, with the curves as in Fig. |^. The lower panel is the same 
comparison for the 1382 jttm, 216 GHz Planck channel. The solid, short-dashed, and dotted 
curves are as before, while the additional long-dashed curve is the simple extrapolation by 
Scott and White (1998). 



-19- 



and 1382 jum, which correspond to the Planck high-frequency channels above and below the 
850/im channel we have considered so far. Again our results indicate that Poisson noise is 
dominant over structure induced by lensing for realistic models of clustering. Note, however, 
that the level of submillimeter anisotropies relative to microwave background anisotropies 
are drastically effected by changes in wavelength. At 545 GHz we expect background 
measurements to be dominated by point source contributions, that is, the far-infrared 
background. At 217 GHz, however, point source contributions are much reduced, becoming 
comparable to microwave anisotropies only at scales smaller than i ~ 20/6> = 1000. As a 
side note we point out that our semi-empirical model shows a much sharper decrease in 
Poisson noise in the 1 to 10 mJy range at 217 GHz than the simple extrapolation by Scott & 
White (1998), due to the sharper turnover of ^ in the TSB model. Thus high-£ microwave 
anisotropy measurements at this frequency may gain more from a careful subtraction of 
point sources than previously expected. 



5. Conclusions 

While the surface density of submillimeter objects was previously unknown to within 
a factor of a thousand, recent observations have given us our first glimpse into the 
submillimeter universe. This breakthrough has allowed authors to quantify the level of 
background emission, confusion limits on future observations, and the effects of lensing on 
the number counts of resolved objects. 

In this work, we have examined the effect of lensing by clusters on the angular 
distribution of submillimeter sources. While lensing does not affect the total flux, it can 
cause a flux decrement of unresolved lensed sources in flux-limited samples. In the case 
of lensing by clusters, this flux decrement is correlated on scales corresponding to the 
cluster correlation length ~ 1.1°, and thus represents a possible contribution to background 
anisotropies at microwave and far-infrared wavelengths. 

Using a simple model for lensing by clusters and a semi-empirical model for the 
evolution of lensed objects, we have quantified this contribution, relating it to the anisotropy 
due to the discreteness of the unresolved source distribution. While slowly evolving cluster 
models predict structure comparable to Poisson noise on angles ~ 1.1°, more likely cluster 
models predict anisotropies to be dominated by Poisson noise. The degree of structure 
induced in an < 1 universe is expected to be more than for the "best guess", f^M = 1 
model, but still fall short of the slowly evolving case we have considered. 

Throughout this work we have assumed a uniform distribution of submillimeter objects 
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and have not addressed the question of the clustering of the sources themselves. This 
clustering length is much shorter than that of galaxy clusters and the clustering weaker, but 
it nevertheless may represent a significant contribution at small angles. At larger angles this 
effect is negligible, however, and thus the issue of structure in the far-infrared background 
itself does not effect the conclusions drawn here. It is with some confidence then that we 
can say that experiments measuring microwave background anisotropies and far-infrared 
background anisotropies due to the clustering of unresolved sources will be able to remove 
all resolved point sources without inducing structure in unresolved sources from lensing by 
galaxy clusters. 



Acknowledgements 

We appreciate valuable input from Rychard Bouwens, Tom Broadhurst, R. Benton 
Metcalf, D. M. Sherfesee, and Edward L. Wright. 



REFERENCES 

Abadi, M. G., Lambas, D. G., & Muriel, H. 1998, ApJ, 507, 526 
Altieri et al. 1998, A&A, accepted, Qastro-ph/9810480D 



Aussel, H., Cesarsky, C. J., Elbaz, D., & Starck, J. L. 1998, A&A, 342, 313 
Balland, C, Silk, J., & Schaeffer, R. 1998, ApJ, 497, 541 
Bahcall, N. A., Fan. X., & Cen, R. 1997, ApJ, 485, L53 
Bahcall, N. A., & Cen R., 1992, ApJ, 398, L81 

Bardeen J. M., Bond, J. R., Kaiser, N., Szalay, A. S. 1986, ApJ, 304, 15 
Barger, A. J. 1998, Nature, 394, 248 

Bersanelli, M. et al. 1996, COBRAS/SAMBA, SCI(96)3, ESA, Paris 
Biller, S. D. et al. 1998, PRL, 80, 2992 
Blain, A. W. 1996, MNRAS 283, 1340 
Blain, A. W. 1997, MNRAS 290, 553 
Blain, A. W., 1998, MNRAS, 297, 511 

Blain, A.W., Ivison, R. J., & Smail, I. 1998, MNRAS, 296, L29 

Blain, A.W., Ivison, R. J., Smail, I., & Kneib J. P. 1998 to appear in "Wide-field surveys in cosmology", 
proc. XIV IAP meeting ( |astro-ph/9806063j ) 

Blain, A. W., & Longair, M. S. 1993, MNRAS, 264, 509 

Blain, A. W., & Longair, M. S. 1996, MNRAS, 279, 847 



-21 - 



Blanchard, A. & Bartlett, J. G. 1998, 332, L49 

Broadhurst, T., Taylor, A. N., & Peacock, J. A. 1995, ApJ, 428, 49 

Carlberg, R. G., Morris, S. L., Yee, H. K. C, & Ellingson, E. 1997, ApJ, 479, L19 

Carroll, S. M., & Press, W. H. 1992, ARAA, 30, 499 

Close, L. M., Hall, P. B., Liu, M. C, Hege, E. K., 1995, ApJ, 452, L9 

Dalton, G. B., Maddox, S. J., Sutherland, W. J., & Efstathiou, G. 1994, MNRAS, 271, L47 

Dyer, C. C, & Roeder R. C. 1973, ApJ, 180, L31 

Eales et al. 1998, ApJ, submitted, ( |astro-ph/9808040| ) 

Eke, V. R., Cole, S., & Frenck, C. S. 1996, MNRAS, 282, 263 

Fan, X., Bahcall, N. A., & Cen, R. 1997, ApJ, 490, L123 (FBC) 

Fixsen, D.J., Dwek, E., Mather, J.C., Bennett, C.L., & Shafer, R.A. 1998, ApJ, 508, 123 

Franceschini, A., Toffolatti, H., Mazzei, P., Danese, L., & de Zotti, G. 1991, A&A, 89, 285 

Franceschini, A., Andreani, P., & Danese L. 1998, MNRAS, 296, 709 

Hauser, M. G. et al. 1998, ApJ, 508, 25 

Henry, J. P., & Arnaud, K. A. 1991, ApJ, 372, 410 

Henry, J. P., Gioia, I. M., Maccacaro, T., Morris, S. L., Stocke, J. T., & Wolter, A. 1992, ApJ, 386, 408 

Holland, W. S. et al. 1999, MNRAS, 303, 659 

Hughes, D. et al. 1998, Nature, 394, 241 

Kcnnicutt, R. C. 1998, ApJ, 498, 541 

Malkan, M. A., & Stecker, F.W. 1998, ApJ, 496, 13 

Navarro, J. F., Frenck, C. S., & White, S. D. M. 1996, ApJ, 462, 563 

Nichol, R. C, Collins, C. A., Guzzo, L., & Lumsden, S. L. 1992, MNRAS, 254, 692 

Peacock, J. A. 1982, MNRAS, 199, 987 

Peacock, J. A. 1986, MNRAS, 223, 113 

Pei, Y. C. 1993a, ApJ, 403, 7 

Pei, Y. C. 1993b, ApJ, 404, 436 

Pei, Y. C. 1995, ApJ, 440, 485 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Puget, J. L. et al. 1996, A&A, 308, L5 

Puget, J. L. et al. 1999, A&A, 345, 29 

Roche, P. F., & Chandler, C. J. 1993, MNRAS, 265, 486 

Sadat, R., Blanchard, A., & Oukbir, J., 1998, A&A, 329, 21 

Sanders, D. B., & Mirabel, I. F. 1996, ARAA, 34, 749 



-22 - 



Saunders, W., Rowan-Robinson, M., Lawrence, A., Efstathiou, G., Kaiser, N, Ellis R. S., & Frenk C. S. 
1990, MNRAS, 242, 318 

Schlegel, D., Finkbeiner, D., & Davis, M. 1998, ApJ, 500, 525 

Scott, D., & White, M. 1998, A&A, accepted, ( |astro-ph/9808003| ) 

Smail, I., Ivison, R. J., & Blain, A. W. 1997, ApJ, 409, L5 

Smail, I., Ivison, R. J., Blain, A. W., & Kneib, J. P. 1998, Invited review at Maryland Astrophysics 
Conference "After the dark ages: when galaxies were young (the Universe at 2 < z < 5)," 
( |astro-ph/981028~ll ) 

Tadros, H., Efstathiou G., & Dalton, G. 1998, MNRAS, 296, 995 

Tan, J. C, Silk, J., & Balland, C. 1999, ApJ, accepted, ( |astro-ph/9904004| ) (TSB) 

Tyson, J. A., & Fischer, P. 1995, ApJ, 466, L55 

Vietri, M., & Ostriker, J. P. 1982, ApJ, 267, 488 

Weinberg, S. 1976, ApJ, 208, LI 

Yee, H. K. C, Ellingson, E., & Carlberg, R. G. 1996, ApJS, 102, 269 

Young, P., Gunn, J., Kristian, J. Oke, J. B., & Westphal, J. A. 1980, ApJ, 241, 507 

Zabludoff, A. L, Geller, M. J., Huchra, J. P., & Ramella, M. 1993, AJ, 106, 1301 



This preprint was prepared with the AAS lATgX macros v4.0. 



